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Abstract. We determined the spatial orientation of the Algol AB close pair 
orbital plane using optical intcrferometry with the CHARA Array, and radio 
interferometry with the European VLBI Network (EVN) . We found the longitude 
of the line of nodes for the close pair being Jli = 48° ± 2° and the mutual 
inclination of the orbital planes of the close and the wide pairs being 95° ± 3°. 
This latter value differs by 5° from the formerly known 100° which would imply 
a very fast inclination variation of the system, not supported by the photometric 
observations. We also investigated the dynamics of the system with numerical 
integration of the equations of motions using our result as an initial condition. 
We found large variations in the inclination of the close pair (its amplitude 
~ 170°) with a period of about 20 millenia. This result is in good agreement with 
the photometrically observed change of amplitude in Algol's primary minimum. 



1 Introduction 

Algol consists of a semi-detached eclipsing binary with an orbital period of 
2.87 days (B8V + K2IV) with an FlIV spectral type star revolving around 
the binary every 680 days. The third component was succesfully observed 
by sp eckle interferometry, and its orbit was precisely determined by iBonneaul 



(11973). This result was refiried b y using the Mark III optical stellar interferom- 



eter (iPan. Shao. Colavital [l99al . 



In the radio regime. [Lestrade et al.l ( 19931 ) detected positional displacement 



during the orbital revolution of the AB pair using the VLBI technique, and 
identi fied the K-subgiant as the sour c e of r adio emission. From their measure- 
ments lKiseleva. Eggleton. Mikkolal (jl998l ) calculated Zm = 100° for the mutual 



inclination of the close and wide orbital planes. This value has been widely ac- 
cepted since then. However, this mutual inclination value cannot be correct 
due to dynamical considerations, because it would produce a fast variation in 
the observable inclination of the eclipsing subsystem, resulting in a fast eclipse 
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depth variati on which contrad icts to the more than century-long photometric 
observations ( SoderhjelmlllQSOl ). 

The aim of this study was to constrain the mutual inclination of the system 
better, requiring the measurement of the longitude of the node for the close 
pair. The other orbital elements are well-known from spectroscopic or photo- 
metric data, but there is a controversy in the value of ^li. Because the expected 
apparent size of the close binary semi-major axis is of the order of 2 milliarcsec- 
onds (mas), we carried out both optical and radio interferometry measurements. 



2 Observations and results 

We carried out interferometric observations both in the optical/near- infrared, 
and in the radio regime. The optical interferometry observations were done 
by the CHARA Array on three nights (2, 3 and 4 December, 2006) in the Kg 
band. The radio measurements were carried out with a subset of the European 
VLBI Network (EVN) on 14-15 December 2006 at 5 GHz. At this observation 
the e-VLBI technique was applied, where the telescopes stream the data to the 
central data processor (JIVE, Dwingeloo, the Netherlands) instead of recording. 
The participating telescopes were Cambridge and Jodrell Bank (UK), Medicina 
(Italy), Onsala (Sweden), Toruh (Poland) and the Westerbork phased array 
(the Netherlands). The measurements and the data reductions were carried out 



accor ding to the standard procedures. (Details can be read in ICsizmadia et al 
l2009l .) 

For the analysis of the CHARA visibilities we developed a model which was 
very close to the one of I Wilson &: Devinnev ( 197ll ) which is based on the Roche- 



model, but it was implemented in IDL to restore the surface intensities into a 
matrix and to calculate the sky-projected picture of the system. Since we had 
only about two dozen visibility measurements, we wanted to limit the number 
of free parameters, choosing just three: the angular size of the semi-major axis 
ai, the surface brightness ratio in Kg band Jrs^ t he angle 17 1 . All ot her 



pa rame t ers we re fixed according to the values given in IWilson et all (|l972l ^ or 



m 



Kim ( 19891 ). After an extended grid search we got oi = 2.28 it 0.02mas 



Jk, = 0.330 ± 0.01, and = 48° ± 2°, respectively. 

The VLBI data were gained during a secondary minimum, which, according 
to our simultane ous optical photometry, occurred at Iq = 2 454084.360 ± 0.003 
( Biro et al. 20071 ). The measurements lasted 9 hours, i.e. ~ 13% of the orbital 



period, but as the projected movement of the source is maximal around the 
minima, the covered projected orbital arc is significantly larger. Due to the 
large inclination, this projected arc is almost a straight line parallel to the nodal 
line, which makes the determination of the longitude of the node more easier. 
Unfortunately, during our measurements a radio flare occured which signiflcantly 
reduced the accuracy of our position determination. Our best fit gave = 
52° lb 3° , but we consider this result less reliable because it is difficult to assess 
how much the radio flare as well as atmospheric phase fluctuations might have 
affected our data (e.g. the displacement of the source during our observation 
was almost twice of the calculated from Keplerian revolution). 

Due to the above mentioned facts, here we concentrate mainly on the 
CHARA results. The true size of the obtained semi-major axis and surface 
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Figure 1. The variation of the inchnation of Algol AB with time. Solid 
line represents the inclinitaion of the close pair observable from the Earth. 
Between the dashed lines the system shows partial eclipses, moreover between 
the dashed-dotted lines the system shows total eclipses for a short time. 



brightness ratio are in very good agreement with previous results, which makes 
our nodal result also plausible. According to the CHARA results, the longitude 
of the node is = 48° ± 2° with an ambiguity of 180°. With the VLBI mea- 
surements we can resolve the ±180° ambiguity and conclude that 0,i = 48° ±2°. 
This is in excellent agreen ient with th e value determined from polarimetric mea- 
surements (r^i = 47° ifc 7° ■ iRudv 19791 ). indicating that polarimetry is an efficient 
tool to determine the spatial orientation of the orbits. 

At this point we can determine the mutual inclination im, and we get = 
95° ± 3°. This value is, however, closer to the exact p er p end i cular i ty than the 
100° which was based on the measurements of Lestrade et al.l ( 19931 ). 



3 Dynamics of the system 



In order to investigate the dynamical behaviour of Algol in the near past and fu- 
ture, we carried out numerical integration of the orbits for the triple system. De- 



tailed description of our code can be found in lBorkovits. Forgacs-Dajka. Sz Regalv 



( 20041 ) . This code simultaneously integrates the equations of the orbital motions 
and the Eulerian equations of stellar rotation. The code also includes stellar dis- 
sipation, but the short time interval of the data allow s that term to be ignored. 
Our input parameters for Algol AB were taken from Kim ( 19891 ) with the ex- 
ception of fli which was set to 48° in accordance with our CHARA result. The 
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orbital elements of the wide orbit were almost identical to the elements given by 
Pan et al. ( 19931 ). As further input parameters, the k2, ks internal structure con- 



stants for the binary members were taken from the tables of iClaret fc Gimene^ 
(119921 '). 



Here we consider only the observable inclination of the eclipsing pair. The 
computed variation of inclination between 7500 BC and 22 500 AD was plotted 
in Fig. [TJ Note that Algol AB does not show eclipses when the inclination is 
lower than 63° or higher than 117°. It shows partial eclipses if the inclination 
is between 63°-117° and moreover, it shows total eclipses when the inclination 
87.3°-92.7. As one can see, Algol started to show eclipses shortly after the 
beginning of the Christian calendar. Total eclipses occured during an almost 
300-year-long interval at the end of the medieval, and at the dawn of the mod- 
ern, scientific e ra. According to our integrations in the time of the discovery as 
a variable star (lMontanarilll67ll ). the inchnation of the close pair was about 



which yielded a total eclipse (with an amplitude of 2.8 mag, when the bright pri- 
mary component is totally eclipsed, and the dominant light source is the distant 
C component), making discovery easier. However, it should be emphasized that 
this time-data are rough approximations only. Since we could not determine the 
position of the node better than ±2°, and for exact calculations one needs an 
accuracy better by one order of magnitude, these numbers should be refined in 
the future. 

Considering the scientific era, one can see that our result suggests an incli- 
nation variation of A^ ~ ~1 ° 6 in t he last century which is in accordance with 
the statement of ISoderhjelml ( 1980l ). 
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